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Abstract: On directed Barabasi-Albert networks with two and seven neigh- 
bours selected by each added site, the Ising model with spin 5 = 1/2 was seen not 
to show a spontaneous magnetisation. Instead, the decay time for flipping of the 
magnetisation followed an Arrhenius law for Metropolis and Glauber algorithms, 
but for Wolff cluster flipping the magnetisation decayed exponentially with time. 
However, on these networks the Ising model spin 5 = 1 was seen to show a sponta- 
neous magnetisation. In this model with spin 5 = 1 a first-order phase transition 
for values of connectivity z = 2 and z = 7 is well defined. On these same networks 
the Potts model with q = 3 and 8 states is now studied through Monte Carlo 
simulations. We have obtained also for q = 3 and 8 states a first-order phase tran- 
sition for values of connectivity z = 2 and z = 7 of the directed Barabasi-Albert 
network. Theses results are different from the results obtained for same model on 
two-dimensional lattices, where for q = 3 the phase transition is of second order, 
while for q = 8 the phase transition is first-order. 
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Introduction 

Sumour and Shabat [U 12] investigated Ising models with spin 5=1/2 
on directed Barabasi-Albert networks 3\ with the usual Glauber dynamics. 
No spontaneous magnetisation was found, in contrast to the case of undi- 
rected Barabasi-Albert networks jH 03 E] where a spontaneous magnetisation 
was found lower a critical temperature which increases logarithmically with 
system size. In S=l/2 systems on undirected, scale-free hierarchical-lattice 
small- world networks conventional and algebraic (Berezinskii-Kosterlitz- 
Thouless) ordering, with finite transition temperatures, have been found. 
Lima and Stauffer [Sj simulated directed square, cubic and hypercubic lat- 
tices in two to five dimensions with heat bath dynamics in order to separate 
the network effects form the effects of directedness. They also compared 
different spin flip algorithms, including cluster flips 0, for Ising-Barabasi- 
Albert networks. They found a freezing-in of the magnetisation similar to 



1 



PU |2] , following an Arrhenius law at least in low dimensions. This lack of 
a spontaneous magnetisation (in the usual sense) is consistent with the fact 
that if on a directed lattice a spin Sj influences spin Si, then spin Si in 
turn does not influence Sj, and there may be no well-defined total energy. 
Thus, they show that for the same scale-free networks, different algorithms 
give different results. More recently, Lima ^U] simulated the the Ising model 
for spin S = 1 on directed Barabasi-Albert network and different from the 
Ising model for spin S = 1/2, the order-disorder phase transition of order 
parameter is well defined in this system. He has obtained a first-order phase 
transition for values of connectivity z = 2 and z = 7 of the directed Barabasi- 
Albert network. Now we study the Potts model for q = 3 and 8 states on 
directed Barabasi- Albert network for values of connectivity z = 2 and z = 7. 
and different from the Ising model for spin 5 = 1/2, the order- disorder phase 
transition of order parameter is well defined in this system. We have obtained 
a first-order phase transition for values of connectivity z = 2 and z = 7 of 
the directed Barabasi- Albert network. 



Model and Simulation 

We consider the Potts model with q = 3 and 8 states, on directed 
Barabasi-Albert Networks, defined by a set of spin variables a taking the 
values 1, 2, and 3 for q = 3 and a = 1,...,8 for q = 8 states situated on every 
site of a directed Barabasi-Albert Networks with N sites. . 

The Potts interation energy is given by 

£ = -./££<W fc (i) 

i k 

where fc-sum runs over all nearest neighbors of i. In this network, each 
new site added to the network selects z already existing sites as neighbours 
influencing it; the newly added spin does not influence these neighbours. To 
study the critical behavior of the model we define the variable e = E/N and 
m = (q ■ maxjfrij] — N)/(q — 1), where rii < N denotes the number of spins 
with " orientation " i — 1,2, and 3 and i — 1,2, ...,8 for q = 3 and q — 8, 
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respectively, in one network configuration. ^From the variable energy we 
can compute the average energy and specific heat and energetic fourth-order 
cumulant, 

u{K) — [< E >) av /N, (2) 
C(K) = K 2 N[< e 2 > - < e > 2 } av , (3) 

HK) = [1 - g^ifU (4) 

where K = J/ksT, with J = 1, and is the Boltzmann constant. Similarly, 
we can derive from the magnetisation measurements the average magnetisa- 
tion, the susceptibility, and the magnetic cumulants, 

m(K) = [< \m\ >] av , (5) 

X (K) = KN[< m 2 >-< \m\ > 2 ] avj (6) 

u *™ = [i-£^]-- (^) 

where < ... > stands for a thermodynamics average and [...] av square brackets 
for a averages over the 20 realizations. 

In the order to verify the order of the transition this model, we apply 
finite-size scaling (FSS) ^Tj. Initially we search for the minima of energetic 
fourth-order cumulant in eq. (4). This quantity gives a qualitative as well as 
a quantitative description of the order the transition [T2]. It is known ^1 
that this parameter takes a minimun value B i rnin at the effective transition 
temperature T C (N). One can show ^3] that for a second-order transition 
liniTv^oo (2/3 — -Bi,mm) = 0, even at T c , while at a first-order transition the 
same limit measures the latent heat |e + — e_|: 

jim (2/3 - B hmm ) - (4 _ e2 _ )2 • (8) 

A more quantitative analysis can be carried out throught the FSS of the 
specific heat C max , the susceptibility maxima Xmax and the minima of the 
Binder parameter B ijrnin . If the hypothesis of a first-order phase transition 
is correct, we should then expect, for large systems sizes, an asymptotic FSS 
behavior of the form ^SJ El , 

Cmax = a c + b c N + ... (9) 
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Bi t min = + bB-iN + ... 



(10) 
(11) 



We have performed Monte Carlo simulation on directed Barabasi-Albert 
networks with values of connectivity z = 2 and 7. For a given z, we used 
systems of size N = 250, 500, 1000, 2000, 4000, and 8000 sites. We waited 
10000 Monte Carlo steps (MCS) to make the system reach the steady state, 
and the time averages were estimated from the next 10000 MCS. In our 
simulations, one MCS is accomplished after all the N spins are updated. 
For all sets of parameters, we have generated 20 distinct networks, and have 
simulated 20 independent runs for each distinct network. 

Results and Discussion 

Our simulations, using the HeatBath algorithm, indicate that the model 
displays a first order phase transition. In plot Fig. 1 we show the dependence 
of the magnetisation M on the temperature T for q = 3 and several system 
sizes. Part (a) shows the curves for z = 2 of top to bottom of N = 250 
to 8000, part (b) the same as part (a) for z — 7. In Fig. 2 we show the 
dependence of the Binder parameter B^K) for connectivity z = 2 and 7 
and various systems size. Part (a) shows the curves for z = 2 of bottom 
to top of N — 250 to 8000, part (b) the same as part (a) for z — 7. The 
Binder parameter clearly goes to a value which is different from 2/3. This 
is a sufficient condition to characterize a first-order transition. In Fig. 3 
we plot the Binder parameter Bi versus 1/N for z = 2 (circle) and z = 7 
(square), and severals systems sizes (N = 250,500,1000, 2000, 4000, and 
8000), Part (a). We show the scaling of the Binder parameter minima, and 
again the first order phase transition is verified. The order of the transitions 
can be confirmed by plotting the values of 2/3 — B itrnin again versus 1/N. For 
a second-order transition the curves goes to zero as we increase the system 
size. Here, the quantity 2/3 — B i min approaches a nonvanishing value in the 
limit of small 1/N, for z = 2 as for z = 7, see Fig. 3, part (b). In the Figs. 
4, 5, and 6 we show the same behavior the Figs. 1, 2, and 3 , respectively, 
but for q — 8. As depicted in Figures 7 an 8, our results for scaling of the 
specific heat and susceptibility are consistent with equations (9,10). In Fig. 
7 we show this behavior for z = 2, part (a), and part (b), z = 7 for q = 3. 
The same occurs with Fig. 8, but for q — 8. 

Conclusion 
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In conclusion, we have presented a very simple equilibrium model on 
directed Barabasi- Albert network |2j. Different from the spin 1/2 Ising 
model, in these networks, the spin q = 3 and 8 Potts model presents a the 
first-order phase transition which occurs in model with connectivity z = 2 
and z = 7 here studied. These results disagree with the results for Potts 
model on two-dimensional lattice ^3], where for q < 4 the transitions is the 
order second and the order first for q > 4 . 
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of this paper. I also acknowledge the Brazilian agency FAPEPI (Teresina- 
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q=3, z=2 




temperature 

Figure 1: Plot of the magnetisation M versus K for q = 3 and severals 
systems sizes (N = 250, 500, 1000, 2000, 4000, and 8000 ). In part (a) z = 2 
and part (b) z — 7. 
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Figure 2: Plot of the Binder cumulant Bi{K) versus K for q = 3 and severals 
systems sizes (N = 250, 500, 1000, 2000, 4000, and 8000 ). In part (a) z = 2 
and part (b) z — 7. 
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Figure 3: Plot of the Binder parameter 2/3 — Bi(K) versus l/N for z = 2 
(circle) and z = 7 (square), part (a) and B imin (K) versus l/N, part (b) for 
q = 3. 
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q=8, z=2 
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Figure 4: Plot of the magnetisation M versus K for q = 8 and severals 
systems sizes (N = 250, 500, 1000, 2000, 4000, and 8000 ). In part (a) z = 2 
and part (b) z — 7. 
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Figure 5: Plot of the Binder cumulant Bi(K) versus K for q = 8 and severals 
systems sizes (N = 250, 500, 1000, 2000, 4000, and 8000 ). In part (a) z = 2 
and part (b) z — 7. 
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q=8, z=2(circle), 7 (square) 
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Figure 6: Plot of the Binder parameter 2/3 — Bi(K) versus l/N for z = 2 
(circle) and z = 7 (square), part (a) and B i:Tnin (K) versus l/N, part (b) for 
g = 8. 
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Figure 7: Plot of the specific heat C max (circle) and susceptibility x 
(square) versus N for q = 3. In part (a) z — 2 and part (b) z — 7. 
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q=8, z=2 
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Figure 8: Plot of the specific heat C max (circle) and susceptibility \. 
(square) versus N for q = 8. In part (a) z — 2 and part (b) z — 7 
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